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Chapter 13. Propagators of the real Klein Gordon field 



13.1 Klein Gordon propagator in position space 

In this section we will study the propagation of the real Klein Gordon field 
from a hypothetical point source In position space we can write. 

/ 1 92 ^2 Q2 g2 ^2^2 X 

Or much more compact, using c = 1 and h = 1 

(□ + m2)V' = C{x^) (13.2) 

Wc define C(-^^) = 5{x^) as a point-source in both space and time. It's 
Fourier transform in momentum space a constant over momentum and 
energy. We can write. 

{-pl + pl+pl+pl + rn)ilj = (V + m2)V = C(p^) (13.3) 

The propagator in momentum space can now be written as. 

^ = 2 C =^ T^P = 2 (13-4) 

The propagator for the Klein Gordon equation in momentum space is 
—ll{q^ — m?), where is the shorthand notation for ^E'^ — p"^ — Py — p^. 

This expression has poles (infinities) in momentum space for certain fre- 
quencies and they need careful attention. In position space we do not get 
infinite values from a plane- wave source. Why? where do they come? Ac- 
tually the propagation from these plane-waves does become infinite but 
only after infinite time. 

The reason behind this is that plane waves stretch out from x is minus in- 
finity to plus infinity. Contributions from farther and farther away regions 
keep coming in and, at the pole frequency, they all add up. The result is 
infinite after t is infinte. Any physical process doesn't continue until t is 
infinite nor is infinite in size. 
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We can split up the physical process of Klein Gordon propagator with the 
following series expansion: 

— 1 /I 
q2 _ ^2 y q2 ^,4 qt> q>i 

Which becomes the following operator in configuration space: 

V{t,x) = D-^ - D-^ + - m^D-^ + .... (13.6) 



Where is the inverse d'Alembertian. This operator spreads the wave 
function out on the Hghtcone as if it was a massless field. The second term 
then retransmits it, opposing the original effect, again purely on the light 
cone. The third term is the second retransmission, et-cetera, ad-infinitum. 

All propagators in this series are on the lightcone. The wave function does 
spread within the light cone because of the retransmission, but it does 
never spread outside the light cone, with superluminal speed. 

We will solve the series first for the 1+ld dimensional case. We may rotate 
the E,p plane and the t,x plane by 45 degrees, they have the same Fourier 
transform correspondence, and use new variables defined by. 



E 



Pi 



V2 ' 



P2 



E + p 
V2 ' 



X 



u 



V2 ' 



w 



t + x 



(13.7) 



1 1 _ 1 1 

E - p E + p 2 pi p2 



(13.8) 



Note that E here ^ ^E with c set to 1. A division by a variable in the 
momentum domain amounts to an integration in the position domain.^ 
We can now express the terms of the series (13.34) as a series of repeated 
integrations of the delta function which originates from the definition of 
the propagator. 

'^Accurately, we should do a convolution with the sign function here while the inte- 
gration is a convolution with the Heaviside step-function 6{x). It docs produce the exact 
result for the forward light-cone. Strictly causal propagation will be studied in detail in 
the succeeding sections 
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{D + m'^)T>{t,x) = 5{t,x) (13.9) 

Writing out the series we get: 
D-^S{t,x) = J J S{t,x) ] dudw = \ e{u)e{w) 

= / / e{u)e{w) ] dudw = I e{u)e{w) u w 

D-^6{t,x) = J J III 9{u)9{w) u w ] dudw = | e{u)e{w) \u^\w'^ 
U-^b{t,x) = J J i[| e{u)e{w)lu^lw'^ ] dudw = ^ 9{u)9{w) ^u^^w^ 

(13.10) 

We can express these terms in the more famihar variable s"^ = — x"^. 
(The t used here is ct with c set to 1). 



uw = (f-x'^)/2= (13.11) 

2 2 

9{u)9{w) = 9^—^) = 9{s^/2) (13.12) 



So we can write for the series: 



Vit, X) = - 9{s-/2) 1 - ^ + - ^ + •••• ) (13.13) 



The series between the brackets becomes a Bessel J function of order zero. 
To see this we use the series expansion for the general Bessel J function. 



fc=0 ^ ' 



In our case we need the zero order Bessel J function: 



13.1 Klein Gordon propagator in position space 



5 



The Bessel Jq series is identical to the one in brackets in the expression for 
the propagator (13.13), so we use it to replace the series: 

V{t,x) = le{sy2)Mms) (13.16) 



We can extend this to any dimensional space using the "inter-dimensional 
operator" The Bessel J function of order zero becomes one of first order in 
3+ Id space. 

= («='^) (13.17) 



Where d is the number spatial dimensions (3 in this case). Prom the series 
expansion for Bessel functions (13.14) we can derive an expression to take 
derivatives of Bessel functions. 

^J^:iM^ = -.-^ j,^,(,) (13.18) 

oz 



Which we can use to apply the intcr-dimensional operator (13.36). As the 
final result we now get for the position space propagator on the forward 
light-cone. 

Real field Klein Gordon propagator in position space 



Vit,r) = eit) ( ^ 6 (.V2) - (.V2) Ji(m.) ) (13. 



19) 



2 J.2 2 
S = t — X 



The inter-dimensional operator is proved in appendix .. In the massless 
limit this expression becomes the Green's function of the classical wave 
equation. The argument used in the delta function (s^/2) is the rigorous 
mathematical version. It can be replaced by one of several equivalent 
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versions with equal volume if the exact shape of the delta function is not 
important. 

The rightmost one explicitly shows the 1/r potential field solution. 



13.2 Real field Klein Gordon forward propagator 

A strictly causal propagator should propagate only within the forward 
light-cone. In the previous section we showed that the standard Klein 
Gordon propagator \I{E^ — — rn^) propagates within the light-cone. 
However, the propagator is symmetrical in E and thus also symmetrical in 
time according to the symmetry-properties of the Fourier transform. 

The propagator is Even in time, 'D{—t, r) = r) We have to find the 
Hubert partner of the momentum-space propagator which represents the 
Odd propagator V{—t,r) = — V{t,r), so that we can obtain: 

^{Even + Odd) = Forward propagator. 
I {Even — Odd) = Backward propagator. 

The antisymmetric Odd propagator is simply the symmetric propagator 
multiplied by the sign function. The derivation of the causal, forward in 
time, Klein Gordon propagator follows closely that of the forward photon 
propagator. We start with the propagator which is symmetric in the energy 
E. (and thus also symmetric in time) 

'I^XliE,p) = E2_p2_^2 (13.21) 

Which corresponds to a propagator in position space which is symmetric in 
time. The triangles in the superscript and subscript symbolize the forward 
and backward light-cone. We can separate this propagator in a positive 
and negative pole. 

1 

E'^ _ p2 _ 
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2y^p2 _)_ j7^2 \e — ■sJjP' + VP? E + sjp^ + m? 

We see that both poles propagate both positive and negative energy plane- 
waves even though their peaks are in different halves. Methods which 
attempt to separate positive and negative frequencies by picking either 
one of the two poles are mathematically incorrect. 

The anti-symmetric propagator is the Hilbert partner of the symmetric 
propagator. The two can be transformed into each other by a convolution 
with i/{'KE). 

* VXl{E,p) (13.23) 

* VZ{E,p) (13.24) 

The square root is just a constant in the above operations. The correct 
correspondence with the massless photon propagator is given by the re- 
placement. 

p <^ + _^ ^2 = (13.25) 

We should keep in mind that the p we are using corresponds to the fre- 
quency domain of r. We could have given it a subscript like in p^. If 
we make the above replacement in equation (??) for the anti-symmetric 
photon propagator then we obtain the anti-symmetric (odd) Klein Gordon 
propagator. 

Vtl{E,p) = -'-^(^S{E-L0p)-d{E + iOp)^ (13.26) 

This propagator propagates only on the mass shell. Once we have the 
symmetric (even) and anti-symmetric (odd) propagator we can combine 
them to get the Forward and backward in time propagators. 



(13.22) 



^-'(--) - (^) 



Forward propagator = ^ ( Even propagator -|- Odd propagator ) 
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Backward propagator = 
^a{E,p) = I ( 

So we obtain the: 

Causal, forward in time, real field Klein Gordon propagator 

ViE.p) = (13.29) 

We see that the modification to make the Klein Gordon propagator causal 
only changes the behavior at the poles, as was the case with the forward 
photon propagator. We reorganize the forward propagator on a pole- by- 
pole base to study the behavior at the poles: 

V\e,p) = (13.30) 

The on-the-mass-shell propagation would be ill defined without the addi- 
tion of the delta-functions which have a magnitude infinitely much higher 
as that of the reciprocal poles alone. See the discussion with Rayleigh's 
theorem (??) 

Somewhat symbolically we can write for the amplitude of the propagation 
if we change the frequency slowly from one side of the pole to the other 
side of the pole. 

—GO — ioo^ — ^ -|- GO (13.31) 

The middle value -ioo'^ symbolizes the contribution of the delta function. 
A real particle with a finite lifetime has a frequency spectrum which is 



VXliE,p) + V':liE,p) j (13.27) 
|( Even propagator - Odd propagator ) 
Vll{E,p) - Vtl{E,p) \ (13.28) 
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not infinitely sharp. The spectrum extends at both sides of the pole. 
There would be destructive interference since the amplitude at both side 
is opposite in sign. 

Minimal changes in frequency would move the center frequency to either 
side of the poles and the destructive interference would disappear. The 
addition of the dominating delta function in the middle of the pole makes 
the on-shell propagation well defined. 

Finally, we also obtain the backward propagator. We do so by simply 
subtract the odd propagator instead of adding it. 

Backward in time, real field Klein Gordon propagator 

T>^{E,p) = (13.32) 



13.3 The KG propagator expanded in the mass term 

We expanded the real Klein Gordon propagator in momentum space into 
the mass term when we calculated the Klein Gordon propagator in 
position space. 



1 2 4 

1 m m 

+ .Z.2 ^2^2 + TB^-Z^ + ■■■■ (13.33) 



In this section we will explore this expansion a bit more. We did see that 
the series can be written in configuration space as. 

T>{t,x) = D-^ - m^D-^ + m^D-^ - D"^ + .... (13.34) 



Where □ -"^ is the inverse d'Alcmbertian operator. The first term of these 
series corresponds to the massless propagator. The subsequent terms can 
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be interpreted as the first, second and third retransmission based on the 
vr? term. 

We can express this series in the inverse d'Alembertian operator as a series 
of Green's functions expressed in the parameter s^, {=i^ — We do so 
by applying the inter-dimensional operator (13.36) on the Id series (13.13) 
to obtain the 3d series. The result is: 

- ^ - 

0{s^/2) ( vn} rn^s^ vp^s^ nfis^ \ 

+ ofioi oi ~ osol A\ + (13.35) 



27r V 22 2^1! 2! 2^2! 3! 283! 4! 



The first term is the Green's function of the massless particle, equal to that 
of the photon. The remaining series can be written as a Bessel J function 
of first order, see (13.19) 

Another intresting observation is that the momentum space series (13.33) 
represents the Id, 3d, 5d, .... Fourier transforms of the Id, 3d, 5d solid 
light cones. A solid light-cone has a constant value everywhere inside the 
light-cone and a zero value outside the light-cone. 

Wc can sec why by repeatedly applying the inter-dimcnsional operator to 
the Green function series in to increase the number of spatial dimensions 
by 2 each step. This comes down simply to a differentiation in per step: 

r>d+2(i,r-) = ^ ^Pd(t,r) (13.36) 

Any term in the form of: 6{s^/2) s^" become constant within the 
light-cone 6{s^ /2) after n steps. This further means that the series ex- 
pressed in the intermediate time/momentum domain represents the Fourier 
transforms of the Id, 3d, 5d ... solid spheres. 

From (??) and (??) and substitution p =^ sgn(p) \/p^ + m? we can derive 
the symmetric in time and Klein Gordon propagator in the t,p domain. 
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As well as the anti-symmetric in time Klein Gordon propagators in the t,p 
domain. 



sin ^a/p^ + tj 
\JjP'-\- VP? 



(13.38) 



We can expand these expression as a series in t. 



sm 



^ , ' = y „ (p^ + m^)" t^n+i (13 39) 



Which does not contain any square roots anymore. The next step is to 
expand these in the mass term vr? which gives us. 



pt -1*3^3 +^^5^5 _lfp7 +^tV... 



+ 



■ijtV +iiV -^tV +|rtV- 



+ 



(13.40) 



+ 
+ 



Where the denominators correspond with the binomial triangle. Each of 
these series can be written as a finite length series in sine and cosine terms. 
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V{t,p) 



pt 



1 2 



hi sin(pt) - ^ cos(pi) 



+ w 



7 sin(pt) - ^ cos(pt) - ^ sin(pt) 



(13.41) 

The first term is the Id Fourier transform of the square-pulse function, 
which can be seen as the Id solid sphere. The second is the 3d Fourier 
transform of the solid sphere. The next series is the 5d Fourier transform 
of the 5d solid sphere, et-cetera. 

The terms between the square brackets are the spherical Bessel functions 
of the first kind which are generated by repeatedly differentiating the sync 
function. 

The first few spherical Bessel functions generated in this way are. 

jo(x) = (13.43) 

X 

, , sin X cos X , . 

n{x) = (13.44) 



smx 




X 




sin a; 


cosx 




X 


3 sin a; 


3cosx 


x^ 





h{x) = ^ - = - — (13.45) 

X'^ X^ X 

The generator of the spherical Bessel functions is thus related to our inter 
dimensional operator from which we can derive the radial Fourier transform 
of a function in any d-dimension from the simplest, one dimensional case. 
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13.4 The square root Hamiltonian 

An operator of interest is the square root Hamiltonian. A direct lineariza- 
tion of the expression. 

-^2— ^ = #2^ = (|pV + mV)v (13.47) 



Where the tildes on the H and p defines them as operator. We can use 
her the series. 

VT+^ = (13.48) 

1 + 1^2 _ 1-4 , J_ ~6 _ _5_^8 , _T_^10 _ ^5512 , _33_^14 _ 429 -16 , 
^ 2^ 8^ ^ 128-f^ ^ 256^ 1024-f^ ^ 2048-^^ 32768-f^ T .... 



This series is found to be represented by a finite expression and we can 
express H in the form of the following expression. 

Q 

ih— tp = Hip = v/ p^c^ + m^c^ ijj = 
ot 

= ± Jl + f (4^)"\i. (13.49) 

\ n!(n- 1)! 22"-i \m?c^ J y ^ ' 

Or written out as a differential operator: 

V = Hip = (13.50) 

{°° (2n-2)! / /i^ \" / 92 ^2 92 \ 
1 ~ X] „!(n-l)! 22n-i i^;;^ J [d^ ^df^d^J J 



Finally, including the EM interactions we get: 

d 

zn— ?/; = Hip =eVip (13.51) 



\ ^ (2n-2)! / 1 Y^l^^^A 
n!(n - 1)! 22"-i I m2c4 ^ 
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In the chapter on the propagators of the complex Klein Gordon fields we 
will find a way to calculate this H via different way and find that the 
application of H is represented by a convolution with a Green's function 
containing Bessel functions. 



